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Recently, some studies have revealed that non-Poissonian statistics of human behaviors stem from 
the hierarchical geographical network structure. On this view, we focus on epidemic spreading in the 
hierarchical geographical networks, and study how two distinct contact patterns (i. e., homogeneous 
time delay (HOTD) and heterogeneous time delay (HETD) associated with geographical distance) 
influence the spreading speed and the variability of outbreaks. We find that, compared with HOTD 
and null model, correlations between time delay and network hierarchy in HETD remarkably slow 
down epidemic spreading, and result in a upward cascading multi-modal phenomenon. Proportion- 
ately, the variability of outbreaks in HETD has the lower value, but several comparable peaks for 
a long time, which makes the long-term prediction of epidemic spreading hard. When a seed (i. 
e., the initial infected node) is from the high layers of networks, epidemic spreading is remarkably 
promoted. Interestingly, distinct trends of variabilities in two contact patterns emerge: high-layer 
seeds in HOTD result in the lower variabilities, the case of HETD is opposite. More importantly, 
the variabilities of high-layer seeds in HETD are much greater than that in HOTD, which implies 
the unpredictability of epidemic spreading in hierarchical geographical networks. 

PACS numbers: 89. 75. He, 87.23. Ge, 05.40.-a 



Since the discovery of non-Poissonian statis- 
tics of human behaviors such as human inter- 
action activities and mobility trajectories, more 
and more scientists have been paying attention to 
the role of these patterns in epidemic spreading. 
Most recent research results showed that both 
time and space activity characteristics respec- 
tively have significant impacts on spreading dy- 
namics. However, it is still unclear to us how the 
spatiotemporal characteristics affect the preva- 
lence. Indeed, the time characteristics of human 
activities is closely related to the space charac- 
teristics. Recently, some studies have revealed 
that non-Poissonian statistics of human behav- 
iors stem from the hierarchical geographical net- 
work structure, in which we investigate how the 
scale-free characteristic of human contact activi- 
ties influences epidemic spreading. We find that, 
compared with homogeneous contact pattern and 
null model, correlations between time delay and 
network hierarchy can remarkably slow down epi- 
demic spreading, and result in a upward cascad- 
ing multi-modal phenomenon. More importantly, 
high-layer seeds arouse large variabilities, while 
low-layer seeds result in several comparable peaks 
of variabilities, which makes the prediction of epi- 
demic spreading hard. This work provides us fur- 
ther understanding and new perspective in the 
effect of spatiotemporal characteristics of human 
activities on epidemic spreading. 
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I. INTRODUCTION 

In modern society, the intrinsic mechanism of epi- 
demic spreading is a noticeable issue. To understand 
the spatiotemporal patterns of epidemic spreading, the 
accurate mathematical models of epidemic spreading are 
used as the basic conceptual tools. There are various 
disease models like SIS (susceptible-infected-susceptible) 
and SIR (susceptible- infected-refractory) [l|, Q ■ With the 
booming development of complex network theory Q , epi- 
demic spreading in complex networks has been strongly 
catching scientists' eyes |3-l3lj. Most studies focus on 
the effect of network structures on spreading dynamics, 
including the small world property ^4-8|, the scale- free 
property [§-[l4j], the community structure [HI HB], and 
the hierarchical structure [l7H24| . etc. Besides, both spa- 
tial distance [25H3G ] and contact capacity [3l|, |32[ were 
found to have nontrivial impacts on epidemic spreading. 

Since the discovery of non-Poissonian statistics of hu- 
man behaviors such as human interaction activities [33j 
and human mobility trajectories (33.[3o|. more and more 
scientists have focused on the role of these patterns in 
epidemic spreading. Most recent research results showed 
that both time and space activity characteristics respec- 
tively have significant impact on spreading dynamics. 
On one hand, the non-Poissonian nature of human in- 
teractions results in slow spreading in the long time 
limit 36-38]. An analytical prediction was then pro- 
posed to understand the emergence of the extremely long 
prevalence time in spreading dynamics f39j . Further in- 
vestigation revealed that this phenomenon mainly stems 
from weight-topology correlations and the bursty activity 
patterns of individuals [40j . By defining the dynamical 
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strength of social ties, an interesting phenomenon was ex- 
plained: although bursts hinder spreading at large scales, 
group conversations favor the local probability of prop- 
agation [4l|. On the other hand, the human mobility 
patterns have a significant influence on epidemic spread- 
ing (42Tj47j . As human traveling statistics follow the 
scaling law, the corresponding simulation results showed 
that the occurrence probability of global outbreaks is 
determined by human travel behavior [42[. Consider- 
ing two distinct individual mobility patterns (i. e., dy- 
namical condensation and object traveling), both theo- 
retical analysis and numerical simulations revealed that 
these patterns have an essential influence on epidemic 
spreading in scale-free networks [I^ |H[ . In the study of 
the fundamental spreading patterns of mobile virus out- 
break, Wang et al. found that a bluetooth virus's spread- 
ing is constrained, which offers ample time for develop- 
ing and deploying countermeasures [45| ■ Recent stud- 
ies presented us that human mobility patterns are often 
dominated by specific locations and recurrent flows [48| . 
Balcan et al. [46[ and Brockmann et al. [47j thus stud- 
ied contagion spreading in bi-directional movements, and 
showed that its dynamics is significantly different from 
random diffusive dynamics. 

Although many researchers have studied the effects of 
time and space characteristics of human behaviors on 
prevalence, it is still unclear to us how the spatiotempo- 
ral characteristics affect the prevalence. Indeed, the time 
characteristics of human activities is closely related to 
the space characteristics. Take human contact activities 
for example. The spatial distance between individuals 
will inevitably lead to the time delay of their contact ac- 
tivities. In this point of view, more attention needs to be 
paid to the effect of spatiotemporal characteristics. Re- 
cently, some studies have revealed that non-Poissonian 
statistics of human behaviors stem from the hierarchical 
geographical network structure (4QL [Hoj , in which we in- 
vestigate how the scale-free characteristic of human con- 
tact activities influences epidemic spreading. We find 
that correlations between time delay and network hierar- 
chy can significantly affect the spreading speed and the 
variability of outbreaks, and it is very difficult to ac- 
curately forecast epidemic spreading in the hierarchical 
geographical networks. 

The paper is organized as follows. In Sec. II, we in- 
troduce the hierarchical geographical network model and 
the propagation processes in two contact patterns. In 
Sec. Ill, we investigate the effects of different contact 
patterns on epidemic spreading. Finally, we draw con- 
clusions in Sec. IV. 



been proposed [50(j . In this model, all nodes are orga- 
nized in L layers. Denote K as the number of first-layer 
nodes, and M as the branching number of the current- 
layer nodes. Each of /th-layer nodes is connected to its 
father node, and two Zth layer node are connected if they 
have the same father node. In the two dimensional plane, 
the whole area is divided into K sub-regions, and K 1 st- 
layer nodes are assigned to locate in the center of them. 
Then, each of the K sub-regions is further divided into M 
sub-sub-regions, with the KM 2nd-layer nodes locating 
in the center. Repeating this process until the Lth-layer 
nodes are generated. Remarkably, there are strong cor- 
relations between the geographical distance and the net- 
work hierarchy in this network: the higher layer a node 
locates in, the farther it is from its father node and the 
other nodes in the same layer. In Fig. [TJ a schematic di- 
agram with L = 3, M = 4, K = 1 is shown. dsA = v2/3 
(the distance between the second layer node B and its 
father node A) is twice of doB = v2/6 (the distance be- 
tween the third layer node D and its father node B), and 
dsc = 2/3 in the second layer is twice of doE = 1/3 in 
the third layer. 



When a random walker continuously jump in this 
network, a power-law-like travel displacement distribu- 
tion P{d) ~ d~ 2 5 is spontaneously generated in the ther- 
modynamic limit t — > oo, where d is the geometric dis- 
tance of a random walker jumping at each time step [50| . 



d=1/6 




II. MODEL INTRODUCTION 

A. Hierarchical Geographical Network Model 

In order to reproduce the scaling law in human tra- 
jectories, a hierarchical geographical network model has 



FIG. 1. (Color online) Illustration of the 2D hierarchical 
structure (L=3, M=4 and K=l) where "diamond", "squares", 
and "solid circles" represent the nodes in the 1st, 2nd and 
3rd layer, respectively. The bolder line of each node repre- 
sents links to its father node and the same layer node on its 
diagonal. Note that "hollow circles" denote no node is in 
these locations. 



3 



B. Time delay in Contact Process 

In various transportation networks, such as railway 
networks and airline networks, the time delay of human 
contact is determined by the geometric distance between 
two neighbors, where time delay is defined as the time 
interval when message is transmitted from one to the 
other. However, in communication networks, such as In- 
ternet and telephone network, time delay of contact has 
no relation with the geometric distance. Thus, contact 
patterns of human activities can be divided into two cat- 
egories. The first pattern is the heterogeneous time delay 
(HETD) of contact activities, which depends on the geo- 
metric distance. For example, cities far from disease ori- 
gins are infected later than those near the origins. Thus, 
the time delay of a contact ry is proportional to the ge- 
ometric distance dij between two neighbors. Obviously, 
the time delay distribution of contact follows power law, 
which is consistent with the distribution of time inter- 
vals between two successive messages arriving at a given 
receiver (5lj . The second pattern is the homogeneous 
time delay (HOTD), in which the distance is ignored. 
For example, in communication networks, computer virus 
spreading is not restrained by the geometric distance. For 
this reason, we set the time delay of a contact a con- 
stant. To compare effects of different contact patterns 
on epidemic spreading, we set the mean time delay of all 
contacts as a fixed value in simulations. It is noted that 
time delay Ty- is generally not an integer such as 1.3. We 
separate its integral and decimal parts as 1.3 = 1 + 0.3, 
and then reset the time delay ry = 1 with probability 
1 — 0.3, while r,j = 1 + 1 with probability 0.3. What 
can be imagined is that these two patterns can result in 
distinct spatiotemporal patterns of epidemic spreading. 



C. Propagation Process 

We study SI (susceptible-infected) [l[ spreading dy- 
namics in contact process (CP) through numerical simu- 
lations. In SI model, '5" and T represents respectively 
the susceptible (healthy) state and the infected state. At 
each time step of contact process, each infected node ran- 
domly contacts one of its neighbors, and then the con- 
tacted neighboring node will be infected with probability 
A if it is in the healthy state, or else it will retain its 
state. To eliminate the stochastic effect of the disease 
transmission, we set A = 1. In simulations, the propa- 
gation processes are as follows: (i) Select a node as the 
initial infected (i. e., seed) and all other nodes are in S 
state, (ii) At each time step, the infected node i in the ac- 
tive state randomly select one of its susceptible neighbors 
j, and then contact node j after time delay ry. During 
this period, node i is inactive, (iii) After Tij, node j is 
infected by node i. Meanwhile, node i is reactivated, (iv) 
The propagation processes will continue until all nodes 
are in / state. 

Based on the above process, we investigate how these 



two patterns influence spreading speed and variability of 
prevalence in CP in hierarchical geographical networks. 
The prevalence is defined as the density of infected in- 
dividuals i{t) at time step t, and the spreading speed 
is defined as new case rate n(t) at time step t, that is 
n(t) — i(t) — i(t — 1). With the spread of epidemic, new 
case rate increases to the maximal value n max , which de- 
notes the occurence of outbreak. In order to analyze the 
impact of the underlying network topology on the pre- 
dictability of epidemic spreading, the variability of out- 
breaks is defined as the relative variation of the preva- 
lence given by [13] 



A [i(t)} = 



m 2 ) - w)y 



(<(*)) 



(i) 



A[i(f)] = denotes all independent dynamics realiza- 
tions are essentially the same, and the prevalence in the 
network is deterministic. Larger A [i(t)] means worse pre- 
dictability that a particular realization is far from average 
over all independent realizations. 



III. THE EFFECT OF DIFFERENT CONTACT 
PATTERNS 

A. Random seed 

In this paper, we study how the different contact 
patterns (i. e., HOTD and HETD) influence epidemic 
spreading. We first pay attention to the case of ran- 
dom seed, that is a node selected randomly as the ini- 
tial seed. Obviously, in Figs. [2] (a) and (b), new cases 
in HETD has the lower peak value n^fJ D « 0.006 < 



,HOTD 



max 

n max " ~ 0.022, and the longer full prevalence time 

rpHETD _ g 00 > rpHOTD ~ 1?5j where the full preva _ 

lence time Tf is defined as the amount of time that all 
nodes of the network are infected by a seed. Meanwhile, 
a upward cascading multi-modal phenomenon in HETD 
is very intriguing, which is consistent with the results 
in Ref . [2J, |45|, l52j . We suppose that this phenomenon 
might origin from the correlations between time delay 
and network hierarchy: time delays of contacts between 
the high-layer nodes are much greater due to longer dis- 
tance from each other. To gain insight into the effect 
of these correlations, we employ null model where time 
delays are randomly exchanged between randomly cho- 
sen links (RNTD), and the time delay- network hierar- 
chy correlations are thus destroyed. In Fig. [5] (a), n(t) 
in RNTD displays almost the same unimodal pattern 
to that in HOTD, except for a slightly lower spread- 
ing (i. e., n™ « 0.020 < n%£ D « 0.022). By 
comparison, we can conclude: although the heterogenous 
time delay of contact activities can result in slow spread- 
ing [H, H3, HU , the very low spreading is dominated by 
the time delay-network hierarchy correlations. More im- 
portantly, it indicates that these correlations certainly 
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lead to the upward cascading multi-modal phenomenon. 
In this network, the outbreaks in low-layer sub-regions 
occur much earlier than that in high-layer sub-regions 
because low-layer nodes are much closer to each other. 
Even if the high-layer nodes are infected, outbreaks will 
not occur in a wider range until their child nodes are 
infected after a long time. Therefore, the peak values 
at t « 7, 33, 78, 170, 350 in Fig. (b) correspond to the 
outbreaks in the 5th, ^th, 3rd, 2nd, and 1st layer sub- 
regions, respectively. 

After that, we also investigate the variabilities of 
outbreaks in the different contact patterns. From 
Figs. [5] (c) and (d), we catch two essential differences 
between HETD and HOTD/RNTD: although the peak 
value of variability Ai max ps 0.45 in HETD is one half 
of Ai max w 0.90 in HOTD/RNTD, there are four com- 
parable peaks in HETD. As pointed out in Ref. [30| . 
the variability changes coincide essentially with the evo- 
lution of the diversity of infected nodes. That is to 
say, the more homogeneously infected nodes distribute 
among all layers, the higher variability is. For the 
case of HOTD/RNTD, the variability is maximal when 
the diversity of layers of infected nodes is the largest 
at t ps 50/60. In addition, there is a small fluctuation 
of the variability at t ps 7 because of the rapid spread 
in the bottom- layer sub-regions (i. e., a M complete 
graph). For the case of HETD, the diversity of infected 
nodes does not vary monotonically over time, because 
there is always a maximum value of the diversity when 
outbreaks occur in the different layer sub-regions. As 
outbreaks successively occur from the low-layer to the 
high-layer sub-regions, the variability displays four com- 
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FIG. 2. (Color online) The evolution of both n(t) and A[i(t)] 
in the different contact patterns. n(t) versus t in HOTD 
(" squares" ) /RNTD ("circles") (a), and HETD ("diamonds") 
(b). Correspondingly, A[i(t)] versus t in HOTD/RNTD (c), 
and HETD (d). The parameters are chosen as N = 9330, L — 
5, M = 6, K = 6. The results are averaged over 10 3 indepen- 
dent realizations. 



parable peaks. It means that the epidemic spreading in 
HETD has the variability over a long period of time, 
which brings a huge challenge to disease control. 



B. Seeds from different layers 

To confirm the dynamical centrality of nodes, it is very 
important to study the effect of different seeds on epi- 
demic spreading [Ha . [HJ . In this section, we investigate 
how seeds from different layers influence the spreading 
speed and the variability of outbreaks. In each simu- 
lation, a randomly chosen node in the designated layer 
is set as the seed. From Figs. [3] (a) and (b), n(t) in 
three different contact patterns follow the same rule that 
the high-layer seeds can accelerate outbreaks. As men- 
tioned in Sec. III. A, the only access to outbreaks in 
the wider range is across the high-layer nodes. Owing 
to time delay-network hierarchy correlations, the high- 
layer seeds doubtlessly make the outbreaks occur much 
fiercer, which corresponds to the larger peak value. It im- 
plies that central cities in transportation networks must 
be crucial regions for pandemic prevention and control. 

Moreover, n(t) in HETD displays an essential differ- 
ence from the other two cases: the peak number of new 
cases rate is equal to the layer I of the initial seed. For 
instance, there is single peak when a 1st layer node is 
set as the seed. On the one hand, owing to the charac- 
teristic of hierarchy, epidemic spreads from the high-layer 
nodes to the low-layer nodes, and thus the infected nodes 
gradually increase with time. On the other hand, the 
early prevalence is too slow due to time delay-network 
hierarchy correlations. After virus outbreaks in the the 
bottom layer, new cases rate rapidly decreases because 
of the effect of network size. Therefore, the seed from 
the I = 1 layer results in only one peak. For the case 
of seeds from the / = 2 layer, epidemic spreading from 
the seed to bottom layer leads to the first outbreak in 
the seed-centered sub-region, and the first peak value is 
greater than the new cases rate value for I = 1 (i. e., 
n; =2 ~ 0.0021 > m = i ~ 0.0018) at t = 120. Owing to 
the farther distance of the seed to both its father node 
and the same-layer nodes, the second peak corresponds 
to the outbreak in the wider range, but the second peak 
value n l ~l x s» 0.0056 < n l =^ x ps 0.0068. Similarly, new 
cases surely have I peaks when the seed is in the Ith 
layer. These results have intriguing implications in hier- 
archical geographical networks: although the high-layer 
seeds make outbreaks occur much fiercer (i. e., the higher 
peak), the low-layer seeds make outbreaks occur much 
earlier (i. e., the earlier peak). 

In Figs. [3] (c) and (d), it is a surprise that the out- 
breaks with the high-layer seeds in HOTD/RNTD have 
the lower variabilities, while the case in HETD is just the 
opposite. For the case of HOTD/RNTD, it takes almost 
the same amount of time to spread upward or downward 
due to no relationship between time delay and network 
hierarchy. When the seed is in the high-layer, epidemic 
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can only spread downward, which has the less optional 
pathways. Thus, the high-layer seeds result in the less 
diversity due to the less optional pathways, which cor- 
responds to the lower variability. As this contact pat- 
tern reflects the characteristics of information diffusion 
in communication networks, the above results demon- 
strate that information from the high-layer nodes have 
two obvious advantages: the faster diffusion rate and the 
better predictability. 

For the case of HETD, it is very difficult for the virus 
to spread upward due to the longer distance from their 
farther nodes. On the contrary, it's much easier to spread 
downward, which can arouse the more optional spread- 
ing pathways and the greater diversity. Thus, the out- 
breaks with the high-layer seeds have the higher variabil- 
ities. For example, A l =* x w 1.35 > A l =£ x « 1.00 > 
^max ~ 0-70. Intriguingly, the maximum value of the 
variabilities in HOTD/RNTD A l =% x rj 0.71 is only half of 
A l =} x fa 1.35 in HETD. It indicates that the time delay- 
network hierarchy correlations make the predictability of 
epidemic spreading worse. On the other hand, although 
the outbreaks with the low-layer seeds have the better 
predictability, several comparable peaks of the variabil- 
ities, which corresponds to several outbreaks in the dif- 
ferent ranges, cause a long-term trouble for pandemic 
prevention and control. 

In order to ensure the universality of the above re- 
sults, other parameters are also chosen to simulate this 
process. Figs. H] (b), (d), and (f) show the results for 
L = 6, M = 4, K = 4. As expected, all simulations 
reveal the same conclusion: the trend of variabilities in 



HETD is completely contrary to that in HOTD/RNTD. 
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FIG. 4. (Color online) In two networks with the different 
parameters, the prevalence and its variability as a function of 
the layer of seed I where "squares", "circles", and "triangles" 
denote the cases of HOTD, RNTD, and HETD, respectively. 
For N = 9330, L = 5, M = 6,K = 6, the peak value of 
new cases rate n max (a), the full prevalence time T/ (c), and 
the peak value of variability Ai ma x (e) versus I. For N = 
5460, L = 6, M = 4, K = 4, n max (b), T f (d), and Aw (f) 
versus I. 



IV. CONCLUSIONS AND DISCUSSIONS 
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FIG. 3. (Color online) The evolution of n(t) and A[i(t)\ for 
seeds from different layers where "squares", "circles", "trian- 
gleups", "triangledowns" , and "diamonds" denote the cases 
of seeds from the 1st, 2nd, 3rd, ^th, and 5th layer, respec- 
tively, nit) versus t in HOTD/RNTD (a), and HETD (b). 
A{i{t)] versus t in HOTD/RNTD (c), and HETD (d). The 
results are averaged over 10 3 independent realizations. 



In conclusions, we have studied the effects of contact 
patterns on epidemic spreading in hierarchical geographi- 
cal networks, and come to a clear understanding that dif- 
ferent contact patterns (i. e., HOTD and HETD) can re- 
markably influence the spreading speed and the variabil- 
ity of outbreaks. First, we focus on the case of the ran- 
dom seed, and find that correlations between time delay 
and network hierarchy in HETD make epidemic spread 
much slower than HOTD/RNTD, and induce a upward 
cascading multi-modal phenomenon. Correspondingly, 
the variability in HETD is lower, but several comparable 
peaks make the long-term prediction of epidemic spread- 
ing hard. Second, we investigate the effect of seeds from 
different layers on epidemic spreading. For three con- 
tact patterns, the high-layer seeds make outbreaks occur 
much fiercer, while the low-layer seeds in HETD make 
outbreaks occur much earlier due to the small distance 
between nodes in the low-layer sub-regions. Interestingly, 
three contact patterns display distinct trends of variabil- 
ities with the different layers of seeds. The high-layer 
seeds in HOTD/RNTD result in the lower variabilities, 
the case of HETD is opposite. It is notable that the 
variabilities of the high-layer seeds in HETD are much 
greater than that in HOTD/RNTD, which implies the 
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unpredictability of epidemic spreading in hierarchical ge- 
ographical networks. To make matters worse, the vari- 
abilities of the low-layer seeds have more comparable 
peaks, which means it is difficult to accurately forecast 
epidemic spreading for a long time. This work provides 
us further understanding and new perspective in the ef- 
fect of spatiotemporal characteristics of human activities 
on epidemic spreading. 
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